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Abstract 

An empirical correlation for predicting critical boiling heat flux (burnout condi- 
tion) under conditions of stable flow for wetting fluids is presented. The corre- 
lation includes data for water, nitrogen, hydrogen, and potassium. Flow passage 
geometries considered include circular tubes, annular tubes, and rectangular 
channels. The effect of swirl devices on critical heat flux is discussed and an ac- 
celeration correlation parameter is included. All fluid property terms are eval- 
uated at the fluid saturation temperature associated with the burnout condition. 

INTRODUCTION 

In the design considerations of boiling fluid devices to be used for the generation of 
large amounts of vapor the best efficiency in terms of large hea. fluxes and low 
wall-to-bulk temperature differences can be achieved within the nucleate boiling 
regime. Transition from the nucleate to the film boiling regime, with its ac- 
companying high wall-to-bulk temperature differences and lower heat fluxes, will 
lower the vapor generation efficiency. The heat flux associated with the transition 
from nucleate to film boiling is herein called the critical heat flux. For a flowing 
system this critical heat flux can be defined as the flux immediately before the 
transition from a high heat transfer coefficient to a lower value at some unspecified 
location along the flow passage. The point along the flow passage surface at which 
this transition in heat transfer coefficient occurs is frequently called the burnout 
point because of the possibility of material failure. The distance measured along 
the flow passage axis from the inlet to the burnout location is designated as the 
critical length. 

Critical heat flux values differing by over 100 percent can easily be obtained in a 
particular research system. These differences frequently can be attributed to: 

(1) the introduction of a compressible volume upstream of the heated test section 
[ 1] , (2) dynamic coupling between the feed system and boiler tube [2] , (3) tech- 
niques used in establishing a critical heat flux condition [ 3] , or (4) gas content of 
the liquid [ 1] . These effects always produce low critical heat flux values. Since 
about 1958, however, many of the preceding effects have been recognized and 
minimized by researchers resulting in better burnout data and understanding of 
two-phase heat transfer problems. 

In the present study, conducted at the NASA Lewis Research Center, a two-step 
correlation of critical boiling heat flux has been developed. The first step con- 
sists of the correlation of the critical heat flux in terms of a critical quality pa- 
rameter for saturated conditions at the boiler inlet (H = 0). The second step is 
the correlation of the increased critical heat flux associated with conditions of 
subcooled liquid at the boiler inlet (H >0). The fluid property parameters de- 
veloped correlate fluids having such widely varying properties as hydrogen, nitro- 
gen, potassium, and water. 
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BURNOUT CORRELATION WITH LIQUID (SUBCOOLING) ENTHALPY 

A typical variation of critical heat flux with liquid (subcooling) enthalpy, H, for 
water is shown in Figure 1 for several flow rates at a fixed passage geometry 
and pressure level (circular tube, 1.08 cm diam; pressure, 6.9 MN/m^). It can 
be seen that the slope of the q c versus H curve increases with increasing values 
of flow rate. Also shown on the curves is the liquid enthalpy point at which the 
critical quality is zero; at greater values of H than that noted, subcooled (bulk) 
burnout occurs while at lower H-values burnout occurs with net vapor generation. 
A single curve can be drawn at each flow rate regardless of whether the critical 
heat flux is in the net vapor generation regime or in the subcooled regime [ 4] . 

The present study showed that the critical quality, X 0 , at which burnout occurs can 
be related to a combination of several well-known heat transfer and hydrodynamic 
numbers and ratios of vapor-to-liquid fluid properties. In very general terms the 
correlation can be expressed as: 


X Q = f(We, Re, Pr, B^_ N ,(i v /^ L , Ap/p v , Dg/E^, L/Dh, a R ) 
The critical quality, X 0 , when H = 0, is defined as: 


( 1 ) 


X D = Qq AW = (4q 0 L)/(GD e \) 


( 2 ) 


A study [ 5] of the various nondimensional relations evolved a dimensionless fluid 
property parameter that, with additional ratios of density and viscosity pi’ovided a 
useful tool for correlation of experimental data. This dimensionless fluid property 
parameter consisted of: 


n b = 


Vs Ap )l P v ( ff ) 1,5 ] 


(3) 


and was arbitrarily designated as a boiling number in reference [ 5] . All fluid 
properties in equation (3) are evaluated at the saturation temperature associated 
with the burnout location. Ng consists of fluid property relations identical to 
those expressed by Reynolds number, Weber number, and the ratio of a character- 
istic dimension of the flow passage to a detachment bubble diameter, d, as follows: 


N B = 


We °h /(ReD h )2 


(Dh/d) 


(4) 


The detachment bubble diameter in equation (4) is that size to which a bubble will 
grow before detaching from a surface with a fluid flow of zero velocity. In terms 
of the Bond -Newton number this bubble size is given by: 


d = [ B 0 - N / f (/3)JV^/g(Ap) (5) 

Equation (5) was simplified for the wetting fluids considered herein by assuming 
that the bubble detachment angle, /3 was sufficiently the same for all the fluids 
and did not vary with surface conditions of the flow passages. It was further as- 
sumed that the term BQ_ N /f(/3) was a constant and equal to 1.0. Equation (5) then 
simplifies to: 
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d =Vcr/g(Ap) (6) 

It was also assumed that the fluid velocity in the region of bubble formation was 
sufficiently low to be neglected as a significant factor in the correlation. Local 
velocity gradients could be included by the use of more complex equations and as- 
sumptions; for example, those of reference [ 6] . The ratio of D^/d inequa- 
tion (4) can now be written as: 

Dh/ d =Vg(^P)Dg/a = Bq_ n (7) 

For the case of nonwetting fluids (mercury) or when wetting agents are used in a 
fluid, consideration of the bubble angle |3 must be included in the present corre- 
lation; however, this is beyond the scope of this paper. 

Mass Velocity 

As a first step, the relationship of the mass velocity to the critical heat flux was 
established. It was pointed out by Rohsenow [7] that a number of researchers 
have determined that the critical heat flux varies with about Vg at low flow rates 
and w.th about G at high flow rates. In evaluating the data ([ 3) , [ 8] , [ 9] ) the 
present author concluded that an exponent e (approximately 0.5 to 1.0) for the 
mass velocity parameter can be related to the channel diameter, the detachment 
bubble size, and fluid properties. A simple equation that satisfies these criteria 



In dimensionless parameters the critical quality then varies with mass velocity 
as follows: 


X o = f ( Re d,v Pr v N B) e ( 9 ) 

Flow Passage Length -to -Size Ratio 

Critical quality obtained with a circular tube [4] (tube diam, 1.08 cm; pressure, 
6.9MN/i. ) is shown in Figure 2 as a function of the mass velocity for severa 1 
L/Djj ratios. This plot, on logarithmic scales, shows that the slopes of the curves 
vary witn L/D^ ratio; at small values of L/D^ the slope is near -1.0 whereas at 
large values of L/D^ the slope approaches -0. 5. The L/D^ relationship demon- 
strated uy these curves can be expressed by: 


(X 0 ) a - 


f( Re d, v Pr v^'B) e 
(L/D h ) b 


(10) 



4 


where the exponent b is a complex variable of density ratio, Ng and L/D^. The 
exponent b is given by: 



The exponent a for the critical quality will be discussed in the next section be- 
cause it is also influenced by flow passage shape. 

Flow Passage Shape 

Critical quality is also influenced by the shape of the flow passage as determined 
from experiments in references [10] , [11] , and [12] . The effect of passage shape 
is accounteo for by another variable in the previously noted exponent a in equa- 
tion (10) as well as a multiplying factor for the right side of equation (10). This 
variable in both cases can be expressed by the ratio Dg/D^. When used as a mul- 
tiplier for the right side of equation (10) the ratio of Dg/D^ is raised to an expo- 
nent of 1. 5. 


The exponent for Xo in equation (10) is now written to include the effects of both 
L/'bh and D e /Dh as follows: 


a = 


1 + (B)(C) 


( 12 ) 


where 


B = 


0. 35*10 


-6 


N- 


B 


1 + (0 0015 Ap/p.,) 2 


1 + (0.00475 L/D h ) 2 - 2 5 


(jn> 


and 


. 0. 125 


C = 


1 + (0.825[(D e /D h ) - l]} 8 
Flow Passage Size 


(14) 


A dimensionless paramtvcr F g , was obtained that related the ratio of the effective 
channel height, b e , [13] to tht detachment bubble diameter in order to account for 
the size trend shown by the data. F s is given by: 


Fs = 


1 + l^d/be) 1 * 67 


(15) 


and is used as a multiplier for the right side of equation (10). 
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For circular tubes and rectangular channels with parallel sides heated b e was 
found to be 0.41 6D or one-half the channel height, respectively. For annular 
tubes with one diameter heated and rectangular channels with only one side heated 
b e is the entire channel height above the heated surface. 

Fluid Properties Parameter 

In addition to the fluid properties included in the previously discussed parameters, 
a separate fluid properties term is necessary to correlate critical quality for dif- 
ferent fluids and for the variation of critical quality with pressure level for a spe- 
cific fluid. A parameter Ff provides the necessary correlation and is given by: 


Ff = 


Pr 


1.5 

v 


(M v /M L )°- 3 (N B xl 06) 0 - 85 


(16) 


The Ff parameter is also a multiplier term for the right side of equation (10). 
Effect of Gravity and Flow Acceleration Devices 


The effect of gravity on critical quality and critical heat iiux has been explored 
both toward zero gravity as well as toward many gravities by numerous research- 
ers. In general it has been found that for "real'* gravity changes, the critical heat 
flux varies with the ratio (g x /g)®' 38 . When radial acceleration devices (such as 
twisted tapes, wire coils, etc.) are used to create an increased gravity field near 
the passage walls, large acceleration forces are required to provide significant 
changes in the value of the critical heat flux and critical quality. Conventional 
methods of calculating radial acceleration values for such devices as tapes and 
coils in tubes yields the following equation [14] : 



The K in t his equation is the ratio of the vapor -to -liquid velocity and is expressed 
as / \Jp L /p v . Thus 


2 7T 2 

1 ::g \ 2 = 

271-2 

p 0 G\ 2 

g(PL *V) D 

\p/dJ 

g (P L p v )D 

[p/dJ 


These relationships do not account for slip between the vapor and the liquid, nor do 
they account for the condition of the liquid on the heated flow passage wall; i. e. , 
whether the wall is fully wetted or rivulet flow is occurring. Consequently, if a^ 
is used in place of g x the effect on critical quality will be greatly exaggerated. 
Water and potassium data (unpublished NASA data and [14] , respectively) show 
that with typical radial acceleration devices g x in terms of a^ can be expressed 
as: 


g x = [l +(0.125 a R )°* 667 ]g 


(19) 
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In the present correlation the effect of radial acceleration is given by: 


F g = (g/g x )°’ 25 


r x ^0.25 

[l + (0.125 a R ) 0 - 667 


( 20 ) 


The Fg-parameter is a multiplier for the right side of equation (10) and a function 
of the exponent e. 


Complete Burnout Correlation for H = 0 

The complete correlation for determining the critical quality and heat flu;: is sum- 
marized in the following relation: 


(X Q ) a = f 


(Re dt yPrvN B Fg) e (D e /D h ) 1 ' 5 (Ff)(F s ) 

(L/D h ) b 


= f(Z) 


( 21 ) 


Representative critical quality data for references [ 3) , [4] , [8] , [9] , [ 10] , ( 11] , 
and [12] are plotted in terms of equation (21) in tho following Figures 3 to 5. It 
should be noted that the low-pressure water data a id that for nitrogen, hydrogen, 
and potassium were all obtained at low flow rates so that the exponent e in equa- 
tion (21) was generally near 0.5. The high-pressure water data, on the other hand, 
was obtained at large flow rates; consequently, the exponent e was effectively 1.0. 


The range of De/D^ for the annular tubes and rectangular channels shown in Fig- 
ure 4 was from 1.12 to 4. 35. The physical size of the annular tubes varied from a 
combination of a heated = 0.95 cm with an unseated D 0 = 1. 8 cm to a heated 
Dj = 1.27 cm with an unheated D 0 - 2.52 cm. The rectangular channels had heated 
widths of 2.54 to 5.34 cm and unheated heights of 0.25 to 1.27 cm. 


The solid curve shown in Figures 3 to 5 was obtained from a data fit using the re- 
lationships expressed by equation (21). The solid curve is expressed as: 


(X 0 ) a = f P - 0035 

1^0. 0035 + 0. 1Z + 


(10Zf 


0. 133 


( 22 ) 


from which 
X 


F" 

° 0 . 


0.0035 


0035 + 0. 1Z + (10Z) 


0. 133/a 


and 


/ GD e \\ 0.0035 1 

V 41 / to. 0035 + 0.1Z + (10Z)J 


0.133/a 


(23) 


( 24 ) 
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CRITICAL HEAT FLUX WITH ENTERING LIQUID ENTHALPY 

The critical heat flux for a fixed geometry and pressure level increases with in- 
creasing subcooling (liquid) enthalpy for a given fluid (Fig. 1). The increase in 
critical heat flux over that given by equation (24) depends on the initial value of 
the critical heat flux when H = 0 and the associated flow rate, L/D^ ratio and 
latent heat; i.e. , critical quality functions. Tl»c present study has assumed that 
the total critical heat flux with liquid subcooling at the boiler entrance can be ob- 
tained from the superposition of q 0 and a Aq that accounts for the effects of the 
liquid subcooling. Thus 


q c = qo + Aq 


(25) 


Experimental data from reference [4] were used to determine the necessary re- 
lationships for Aq. A good representation of these data can be obtained from the 
following equation: 


Aq = 


1 5500 ( ^]° ' | til | ' 

\ P v/ UJ 


67, 


\0. 33 




^TO^(HA) 


+ 


8.5(1 - 



1 + (0.002 X 0 L/D h ) 2 (Pv/Ap) 0 - 33 (M L /Mv) 


(26) 

The value of Xq in equation (26) is obtained by use of equation (23). Comparisons 
of the measured Aq for representative flow conditions and circular tubes are 
shown in Figure 6 as a function of subcooling (liquid) enthalpy. Similar compari- 
sons for a representative annular tube (Dq/Dj = 1.8 cm/0.95 cm) and a rectangular 
channel (0. 64*5. 33 cm) are shewn in Figure 7. While the data shown in Figures 6 
and 7 are all for a pressure of 6.9 MN/m 2 , similar results were obtained over a 
range of pressures from 3.85 to 13.8 MN/m 2 . 

CONCLUDING REMARKS 


From an analysis of critical heat flux data obtained with water, nitrogen, hydrogen, 
and potassium, empirical equations have been evolved that correlate critical quality 
and critical heat flux over a wide range of operating conditions and flow passage 
sizes and shapes. Correlation was achieved by using established dimensionless 
groupings and fluid property ratios (vapor to liquid) known to influence boiling. All 
fluid properties are evaluated at the saturation temperature associated with the 
burnout location. The correlation is suitable for both net vapor generation condi- 
tions as well as bulk (subcooled) burnout conditions. The inclusion of an accelera- 
tion term should be useful in predicting the improvement in boiler burnout by the 
application of internal flow swirl devices. 

On the basis of a preliminary study, the author believes that critical heat flux val- 
ues for more complex geometries, such as rod bundles, may be estimated by use 
of the present correlation. Such estimates require the selection of suitable shape 
and size factors. For example, a single rod inside a rod bundle can perhaps be 
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analyzed by considering the flow passage periphery to be represented locally by 
an annulus or channel heated from both sides. A rod on the outside of the bundle 
adjacent to the shell can be considered similarly but with only one surface heated. 

NOMENCLATURE 


ag radial acceleration 

b e effective channel height 

Bq_n Bond Newton number 

B' XT modified Bond -Newton number 

O-N 

D flow passage diameter 

d detachment bubble diameter 

hydraulic diameter (based on 
wetted perimeter) 

D e equivalent diameter (based on 

heated perimeter) 

G mass velocity 

g acceleration at standard earth 

gravity 

g x acceleration at other than 

standard earth gravity 
H liquid (subcooling) enthalpy 

measured from saturation 
L channel length to burnout 

Ng boiling number 

P boiler insert pitch 

Pr Prandtl number 

Q 0 rate of critical heat transfer 

q c total critical heat flux, 

q G + Aq 

q Q critical heat flux with satu- 

rated fluid at boiler inlet 
Aq increase in critical heat flux 

above q G when fluid at 
boiler inlet is subcooled 
(H >0) 


Re Reynolds number 

W mass flow rate 

We Weber number 

X thermodynamic quality 

X Q critical quality with saturated 

fluid at boiler inlet (H - 0) 

Greek Symbols 

0 bubble contact angle 

\ latent heat 

H dynamic viscosity 

p density 

A p density difference between 

liquid and vapor 

a vapor-liquid surface tension 

Subscripts 

d bubble diameter 

hydraulic diameter 
i,o inside and outside diameters of 
annular tube, respectively 
L liquid 

v vapor 
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Figure 1. - Typical rariation of 
critical heat flux with liquid 
(subcooled) enthalpy Ql]. 



GxlO' 3 , kg/sec-m 2 

Figure 2. - Typical variation of 
critical quality with mass velocity 
(H = 0), [4]. 
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Figure 4. - Correlation of critical quality for channels of various cross sections. 
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Figure 5. - Correlation of critical 
quality in circular tubes for vari- 
ous fluids. Pressure, <0.35 
MN/m 2 . 
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GEOMETRY 
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L/D h 

0 

ANNULAR 
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11] 

a 

ANNULAR 

.363 

209 

M 

0 

RECTANGULAR 
EQ. (26) 

.498 

83 
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Figure 6. - Correlation of liquid 
enthalpy effect on critical heat 
flux. Circular tubes; pres- 
sure, G.9 MN/m 2 , £ 4 ]. 


Figure 7. - Correlation of liquid 
enthalpy effect on critical heat 
flux. Various shaped channels; 
pressure, 6.9 MN/m 2 . 
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